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Abstract 

For a right-invariant system on a compact Lie group G, I present two methods to 
design a control to drive the state from the identity to any element of the group. The 
first method, under appropriate assumptions, achieves exact control to the target 
but requires estimation of the 'size' of a neighborhood of the identity in G. The 
second method, does not involve any mathematical difficulty, and obtains control to 
a desired target with arbitrary accuracy. A third method is then given combining 
the main ideas of the previous methods. This is also very simple in its formulation 
and turns out to be generically more efficient as illustrated by one of the examples 
we consider. 

The methods described in the paper provide arbitrary constructive control for 
any right-invariant system on a compact Lie group. I give examples including closed 
multilevel quantum systems and lossless electrical networks. In particular, the re- 
sults can be applied to the coherent control of general multilevel quantum systems. 
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1 The Lie algebra rank condition of geometric con- 
trol theory 

Consider a control system of the form 

x = f{x,u), (1) 

where x is the state varying on a compact Lie group and u the controL The system is said 
to be right invariant if, denoting by x(t,u,s) the solution of (1) corresponding to initial 
condition s and control function u, we have 

x{t,u, s) — x{t,u,l) o s, (2) 

where 1 denotes the identity of the group and o is the multiplication of the group. To 
be concrete, we shall consider the case of matrix groups where the group operation is the 
standard matrix multiplication, with particular attention to subgroups of SU{n), given 
the potential application to quantum systems. In particular, we shall consider systems of 
the form 

X = A{u)X, X(0) = 1, (3) 

where 1 is the identity matrix and the matrix A{u) is in the Lie algebra associated with 
G for every value of the control u. This equation models many systems of interest. In 
particular closed (i.e., not interacting with the environment) finite dimensional quantum 
systems which are coherently controlled (i.e., through a variation of their Hamiltonian) 
are modeled this way. In this case, equation (3) is Schrodinger equation. We refer to [7] 
and references therein for several examples and introductory notions on Lie groups and 
Lie algebras in the context of quantum control. 

If we restrict ourselves to piecewise constant controls, the problem of control for sys- 
tems (3) can be described as follows. Assume that we have a linearly independent set of 
matrices 

T:={A,,...,A^}. (4) 
To each of them there corresponds a semigroup 

Sr.= {e^''\t>0}, i = l,...,m. (5) 

The problem of control to a matrix Xf is to choose N elements Xi, I = 1, . . . ,N G Sj, 
for some j = 1, . . . , m, such that HiLi Xi = Xf. If such elements exist Xf is said to be 
reachable. The question of the set of reachable matrices is a standard one in geometric 
control theory. The result in the following Theorem 1, known as the Lie algebra rank 
condition, is classical [11] and provides the answer for compact Lie groups. 

Let C be the Lie algebra generated by the elements in JF defined as the smallest Lie 
algebra containing and denote by the connected Lie group associated with C. We 
shall call C the dynamical Lie algebra associated to the system. 



2 



Theorem 1 [11] Consider the Lie group and assume it is compact. Then, the set of 
reachable values for X in (3) is equal to e'-' . 

This result has been elaborated upon in several papers and applied to quantum me- 
chanical systems (cf. [1], [7], [10], [18]). In particular, in the case of (closed) quantum 
mechanical systems £ is a subalgebra of the unitary Lie algebra u{n) and, as such, can be 
written as the direct sum of an Abelian subalgebra and a semisimple subalgebra to which 
there corresponds a compact Lie group. That is, modulo an Abelian subgroup which 
commutes with all of e'^, is compact (cf. [6] and [17]). In particular, is compact 
if £ = u{n) ox L — su{n) in which case, the system is called controllable and is the 
group of unitary matrices U (n) or special unitary matrices SU (n) , respectively. 

The original proof given in [11] is not constructive, i.e., in our setting, it does not show 
how to alternate elements in the semigroups Sj in (5) to obtain a given target Xf G e^. 
We show how to obtain this in two ways in the following two sections. The main ideas 
are then combined in a third method in section 4. The first method, described in section 
2, achieves exact control if the subgroups corresponding to the semigroups in (5), i.e., 

Sj := {e^^'\t e R}, j = l,...,m, (6) 

are closed. Otherwise it obtains control with arbitrary accuracy as it follows from Propo- 
sition 2.1 and Remark 2.2 below. This proposition allows us to replace an exponential of 
the form e^* with t < with an exponential of the form e^* with t > which approximates 
it with arbitrary accuracy. This result will be utilized for the following two methods as 
well. 

2 Method 1: Exact constructive controllability 

The method we are going to describe is a consequence of the proof of the Lie algebra rank 
condition. Theorem 1, given in [7] and the result on uniform finite generation of compact 
Lie groups given in [5] . Let Xf e be the target state. We want to show a way to obtain 
Xf as a product of elements in (5), if not exactly, at least, with arbitrary accuracy. We 
are first going to relax the problem by allowing the use of elements in the subgroups (6) 
rather than only elements of the semigroups (5). We shall show later how to overcome 
this problem (see Proposition 2.1 and Remark 2.2). 

Since is compact the exponential map is surjectivc, that is, there exists a matrix 
A E C such that e"^ = Xf, for every Xf} This also implies that, given any neighborhood 

K of the identity in e^, wc can choose an integer M sufficiently large such that cm = 
1 

Xf^ e K. Now, assume first that JF is a basis for C, that is, no Lie bracket is necessary 
to obtain a basis of JC. This implies that, by varying ii, . . . , in a neighborhood of the 
origin in R^, K := {X = e^sts ^Am-itm-i . . . ^Mti |^^^ ^ g gjygg ^ neighborhood 

^See, e.g., [14] and [16] for a study on the generalization of this result. See also [9] (Theorem 6.4.15) 
for the theorem on existence of the logarithm of a matrix. 
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of the identity in e"^ and, in particular, it contains for sufficiently large M. That is, 
we can find real values t\, . . . ,tm such that 



Therefore, by using elements from the subgroups (6) we can obtain e^. Now assume T 
is not a basis for C Since J- :— {Ai, . . . , A^} generates all of C, there exist two values 
1 < k,l < m such that the commutator [Ai, Ak] is linearly independent of {Ai, . . . , Am}- 
This implies that there exists a value t e R such that F := e^'M^e"^'* is also linearly 
independent. To see this, assume it is not true and write e'^'^^lfee"^'* as 

m 

e^'*A,e-^'* = ^a,(i)A„ (8) 

for every t. Taking the derivative with respect to t ai t — 0, gives [Ai, A^] — J^JLi Qi(0)^j, 
which contradicts the fact that [Ai, A^] is linearly independent of {^i, . . . , A^}. Let t be 
such that 

F := e^'*Afce-^'* (9) 

We can add F to {Ai, . . . , A^} and still have a linearly independent set. Moreover, we 
can express every exponential e^* in terms of exponentials of Ai and A^ since e^* — 
^Ait^Akt^-Aii^ Define A^^+i ■— F. If {^41, . . . , A^, ^m+i} is a basis of C then we can pro- 
ceed as above and obtain a neighborhood of the identity in by varying {ti, . . . tm+i} £ 
R"'+^ Such a neighborhood is given by K := {H^^/ e'^^'^lti, . . . ,t„,+i G R }. If that 
is not the case, then we observe that {Ai, . . . ,^^+1} is still a set of generators for C 
and, as above, there must exist two elements Aj. and Ai in {Ai, . . . , A^^i}, such that 
[74fe,74;] is linearly independent of {Ai, . . . , A^+i} and therefore for some t, ^^+2 '■— 
e'^'*Afee~^'* is linearly independent of {^41, . . . , A^+i}. The exponential 6^'"+^* again 
can be expressed in terms of exponentials of Ai, . . . , A„i+i and therefore in terms of 
exponentials of Ai, . . . ,Am- Proceeding this way, one finds dim(£) — m new matrices, 
{Am+i, • • • , ^dim(/;)} which together with {Ai, . . . , A„i} form a basis for C. By 

taking Hj™^^'' e'^-'*^ with tj G R, j = 1, . . . , dim(£), we obtain all the elements in a 
neighborhood of the identity and in particular cm. Repeating the sequence M times we 
obtain e^. 



In the expression of cm and therefore in the expression of e , there will be some 
exponentials with negative t, i.e., some elements in the subgroups (6) which are (possibly) 
not in the semigroups (5). There are ways to minimize the number of these elements in 
the full product, for example by placing together matrices which come from similarity 
transformations with the same matrix so as to have cancelations of the type e'^^^'^e^'^^^^ = 
gAj(ii-t2) Also, in many cases, the orbits {e^^*|t e R } are periodic (closed), which allows 
us to assume all the tj^s positive, without loss of generality. However, if this is not the 
case we can use the following fact. 
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Proposition 2.1 Let e ^1*' an element of a compact Lie group e"^. For every e > there 
exists a,t>0 such that^ 

||e-^l*l - e^*|| < e. (10) 

Proof. Consider e~^l*l and the sequence e'*^!*!, which by compactness of has a con- 
verging subsequence e"'^'^)^'*'. We have hm^^oo e*^"'^^"'"^^""^'^^""^^'^'*' = e~-^l*l. Therefore 
there is k such that ^\eiMk+i-)-n{k)-i)B\t\ _ g-s|t||| ^ ^j^g proposition holds with 

t ^ {n{k + 1) - n{k) - l)\t\. □ 



Remark 2.2 The proof given above follows the one given in [11]. A different, more 
concrete, proof can be given for Lie subgroups of U{n), which is the case that interests 
us the most. In that case, using the Frobenius norm of matrices, we have 



n-Y,cos{uj{t+\t\)), (11) 

where icUj, j = 1, . . . ,n are the eigenvalues (possibly repeated) of B. If we can choose 
i > so that ^ 

[l-cos(a;,-(t + |t|))]<^, (12) 

for every j = 1, . . . ,n, then (10) is certainly satisfied. If g := arccos (l — then, we 
satisfy condition (12) if we are able to find i and integers mj, j = 1, . . . ,n such that 

+ \t\) - 2'Kmj\ < g. (13) 

However, according to Dirichlet's approximation theorem (see, e.g., [2]), given a natural 
number N and n reals ai, . . . , «„, we can find positive integers a,bi, . . . ,bn, with 1 < a < 
N'^ so that \aja — hj\ < ^. This result can be apphed to satisfy condition (13) identifying 

aj with and choosing jf < ^ and choosing nij = bj and t so that = a. Notice 
that since a > 1, t > as desired. For the problem to find a and bj^s, there are several 
algorithms in the literature (cf. [12] and [13]). Notice, in any case, that we are only 
interested in a, which determines i, and since a is bounded from above by A^", it can be 
always found, in principle, by exhaustive search. 



eBf_e-B|t| 



V2 



\ 



We can summarize the given method as follows: 

1. Given JF := {Ai,...,A„i} find, via similarity transformations, dim>C — m more 
matrices {74^+1, . . . , ^dim(£)} so that {^41, . . . , ^dim(£)} is a basis for C 

2. Take the (principal) logarithm of Xf, A, so that = Xf. 

^Whenever we do specific computations involving norms of matrices we use the Frobenius norm || A|| := 
^/Trace{AA^). 



5 



3. Find M (sufficiently large) and h, . . . , tdimic), so that 

dim(/;) 

ThenX; = e^=(n-r/''^e-^^*^)''. 

4. Replace the exponentials of the matrices Am+i, . . . , Adim(£) with expressions involv- 
ing the exponentials of {Ai, . . . , A^} as obtained from step 1. 

5. Replace every exponential c^*, (B G JF) involving negative t with its approximation 
involving positive t. This can be obtained with arbitrary accuracy according to 
Proposition 2.1 and Remark 2.2. 



Remark 2.3 In the above procedure, step 3. is decidedly the most difficult one since it 
requires the solution of nonlinear equations involving the exponentials of matrices. The 
solution is guaranteed to exist for M sufficiently large. This task is obviously easier for 
low dimensional systems. It must be remarked however that there is some flexibility in 
the choice of the matrices A^^i, . . . , Adim(/:)) because of the choice of the pair A^, Ai and 
of the times i (cf. (9)). We can use this flexibility to make these matrices as simple as 
possible (e.g., block diagonal, sparse, etc.) so that calculating the exponential is easier. 
Another type of flexibility, which may be used in calculations, is the fact that the way 
exponentials are arranged in (14) is arbitrary. Any different order will give a neighborhood 
of the identity also. The methods described in the following two sections do not present 
this problem. 

Remark 2.4 The last step of the method can be achieved exactly (i.e., without involving 
an approximation) if the orbit associated with the given matrices JF := {Ai, . . . , A^} are 
periodic. In this respect, notice that, if this is the case, all the other matrices obtained 
by the method also have associated periodic orbits (their eigenvalues are the same as the 
ones of the original matrices). Therefore, for a given matrix B, and negative t, we can 
choose a positive t, such that e-^* = e^*. 

Remark 2.5 [5] It is interesting to give an upper bound to the number of exponentials 
involved in obtaining a neighborhood of the identity according to the described method. 
Let us assume that, at every step, we only produce one new linearly independent matrix. 
For the given matrices {Ai, . . . ,Am}, we need only one exponential, but for the matrix 
obtained at step 1 we need three exponentials. In general, at step j, j > 2, the worst 
case scenario is when we combine a matrix obtained at step j — 1 (giving the similarity 
transformation {Ai in (9)), which requires dj-i exponentials, with a matrix obtained at 
step j — 2, which requires exponentials. The total number of exponentials at step j 
is therefore dj = 2dj_i + dj_2- Therefore having deflned recursively the numbers dj as 

do = 1, = 3, dj = 2dj_i + dj-2: (15) 
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the number of exponentials required is 



dim£— m 



2.1 Example 



We illustrate this method with a simple example of the quantum control of a two level 
system, i.e., a control problem on SU{2), which is compact. Recall the definition of the 
Pauli matrices 

"^'■^{l o)' "y'-^i-i o)' ^^^=(0 -i)- 

Let ^ := {Ai, A2}, with Ai := iaz and A2 := i{ax + ay). Calculate e^^*y42e~^^* which 
for t — —^71 gives A3 — —iy/2ay, which is linearly independent of Ai and A2, and along 
with them it forms a basis of su{2). A straightforward calculation gives 

M_fe''' 0\ A,t,_( cos(v^t2) e^^ sin(v^t2)\ 

"VO e-^*ij' " V-e-^^sin(v/2t2) cos(v^t2) J ^ ' 



- sin( V2t3) cos (72^3) 



and the the set 




Si,2,3 := {e^'''e^'''e^'''\h,t2,ts e R }, (19) 
covers a neighborhood of the identity in SU{2). Assume now our target state Xf is 

(20) 

We first try to see if Xf is in the set 5*1,2.3 in (19). Therefore we must be able to choose 
ti and ts so that P := e~^^^^Xfe~^'-^^^ has the form e^^*^. This means in particular that 
the difference between the phases of the Pi^2 element and Pi^i elements in P is ^. As 

a straightforward calculation shows, Pi,2Pi*i = | independently of the choice of ti and 

1 

^3. Therefore Xf ^ 5'i,2,3. We replace Xf with Xf. The same calculation shows that, 

for every ti, Pi,2Pi*i = ^sm{2\/2ts) + and, therefore, the choice ^3 := ^ achieves 
the desired phase difference. Then, we can choose ti to impose that the element Pi,i has 
phase zero (it is real). This leads to ti — With these choices, we have 

Ci 1 \ 

i^_-3^ ^1 I • (21) 



Comparing this with e^*^ in (18) leads to the choice t2 — With these choices 
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In terms of the original available matrices, Ai and A2, we have 

Xf = (e^i*ie^^*^e-^i^e^^*^e^i^)' , (22) 

where ^1,^2,^3 are the ones found above. The presence of the negative 'time' — ^ in 
the third exponential, does not pose any problems since the one dimensional subgroup 
associated with Ai (as well as any other matrix in su{2)) is periodic. 

A similar treatment shows that, had we chosen to work with the set 

^1,3,2 {e^i^i 6^3*^6^^*^ 1^1, t2,heR}, (23) 

we would have achieved Xf with just three exponentials. This shows that the order in 
which the exponentials are chosen may be important. 

It must be said that for the special case of SU{2) there are many more techniques 
which may be preferable to the one advocated here. For example, since one has available 
both icTz and iay one could have applied a simple Eulcr decomposition. In general it is also 
possible, for general target matrices, to find the factorization with the minimum number 
of factors [4]. Our goal here was to illustrate the method on a simple, easily computable, 
case. We remark that even for large dimensional Lie groups, one can combine these ideas 
with Lie group decompositions for which there exists a large set of tools [7] . 

3 Method 2: Constructive controllability with arbi- 
trarily small error 

In this and the following section we illustrate methods which do not require the solution 
of nonlinear algebraic equations, such as (7), but can be implemented with simple linear 
algebraic techniques. The algorithms achieve control to the target with arbitrary small 
error. 

Reconsider the available set of matrices in (4). As before, we relax the requirement 
to use only elements in the semigroups (5) and use elements in the subgroups (6). We 
can then replace elements in the subgroups with elements in the semigroups as done in 
the previous section. We start with a definition 

Definition 3.1 A matrix H is said to be simulable with the set J-' if there exist r contin- 
uous, strictly increasing, functions fj, j — 1, . . . ,r, with fj{0) — 0, defined in an interval 
[0, e), such that 

^Hx ^ -Q ^Ljfj(x) ^ o(x^+^), (24) 
i=i 

for some matrices Lj e IJ — and^ a 5 > 0. 
^—J^ denotes the set {—Ai,—A2, . . . , —Am}- 
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If a matrix H is simulable, we can control from the identity to with the desired 
accuracy using elements in the subgroups (6) (and therefore of the semigroups (5)). 



Lemma 3.2 Assume (24) holds. Then 



lim n 

77. — >nc> I 



Proof. If (24) holds then 



(25) 



hm TT 



lim 

n— »oo 



n 



1+5 



(26) 



However, we have this standard limit in matrix analysis (see, [9] section 6.5) 

1 



lim 

n— »oo 



O 



1+5 



n 



which proves the lemma. 



(27) 
□ 



From the point of view of constructive controllability, this lemma says that, for each 
simulable H, we can put together a product of exponentials of elements in which, 
repeated a large enough number of times, approximates, with arbitrary accuracy, . 



Theorem 2 Every H in the dynamical Lie algebra C is simulable. 

Remark 3.3 This theorem along with Lemma 3.2 and Proposition 2.1 give an alternative 
proof of a slightly weaker form of the Lie algebra rank condition of Theorem 1. Since e'~ 
is compact, for every Xf in c^, there exists an H E jC such that = Xf. Theorem 2 
and Lemma 3.2 say that we can find a sequence of reachable points converging to Xf for 
every Xf. Therefore the set of reachable states is dense in e^. 

Remark 3.4 Elaborating on the proof of the Theorem 2, wc will also show how to choose 
the elements Lj E J-'[j —T and the functions in (24) so as to make the controllability 
result constructive. We shall discuss this after the proof. 

Proof. The proof is similar to the one given in [3] in the context of quantum walks 
dynamics. In particular, we will show that the set of simulable elements H isa Lie algebra 
containing and this will be sufficient since C is the smallest Lie algebra containing JF, 
by definition. 

First of all, it is clear that every element in JF is simulable, since equation (24) holds 
with r — 1 and = 0. Therefore the set of simulable matrices contains 
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Moreover if H satisfies equation (24), tfien we fiave 

1 1 

g-i/x _ -Q ^-Ljf,{x) _ -Q ^-Ljf,(x)Q(^^l+S^^-Hx^ (28) 
j=r j=r 

and by expanding the exponentials it follows that the last term is also an 0{x^'^^). There- 
fore —H is also simulable. Moreover, for a > 0, (24) holds for aH with fj{x) replaced by 
fj{ax) and 0{x^+^) replaced by 0{a^+^x^+^) = 0{x^+^). If (24) holds for Hi and H2, i.e., 
we have 

^H,x ^ Y[ e^f^^ + Oi(x^+^0, ^ = 1, 2, (29) 

combining this with 

^ Q(^2) ^ e^l=^e^2^ (30) 

gives^ 

^{H,+H,). ^ -Q gL|/|(x.) -Q ^L]f]i.) ^ 0{X^^'), (31) 

with 5 = min{5i, ^2, 1}. Therefore, if Hi and i?2 are simulable, so is Hi + H2. These 
arguments show that the set of simulable matrices is a vector space. 

To show that it is also a Lie algebra, we have to show that if Hi and H2 are both 
simulable so is [Hi, H2]. In order to see that, write (29) in the form 

e^i* = Ti(t) + Oi(f+''), e"^' = T2{t) + 02{t^^''), (32) 

i.e., by replacing the products with the functions Ti and T2. This also gives (cf. (28)) 

e-^i* = T-\t) - Ti\t)Oi{t^+^')e-"'\ e"^^* = ^^^(i) - T2\t)02{t^^^')e-"^K (33) 

We use the exponential formula (sec, e.g., [9] Section 6.5) 

Using (32) and (33) in (34), we have 
^[H,,H,]e ^ (5(^3) ^ i^^-i _ TfiOie-^i*) {T:^^ - T^^02e-''^') {Ti + d) {T2 + O2) ■ (35) 

Expanding the right hand side, omitting terms that are clearly 0{t"), a > 2, since they 
contain the product of two O functions, we have 

^[H^,H2]t^ + 0{t^) = T{^T^^TiT2 + T{^T^^Ti02 + T{^T^^0iT2 (36) 
-Ti^T2^02e~"^*TiT2 + Ti^Oie-"^*T2^TiT2 + 0{e). 



^Here and elsewhere, we use the notation O for a generic O-function and we use indexes Uke in 0\ 
and O2 when we want to highlight a particular O-function. 
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Expanding in McLaurin series the functions multiplying the Oi and O2, we see that the 
terms corresponding to the first terms of the expansion cancel, leaving only terms of the 
form 0{t^ ) with (3' > 2. In conclusion, we have 

^{HuH2]t^ ^ T-\t)T^\t)T,{t)T2{t) + 0{t^), p>2, (37) 

and by setting t = ^/x, we obtain 

^[H,,H,]. ^ T{\^)T-\^)T^{^)T2{V^) + O(x^), /3 > 0, (38) 

which shows that [Hi , H2\ is simulable as well, and completes the proof. □ 

In order to use Lemma 3.2 and Theorem 2 for control, we need to show, given iJ, how 
to find the matrices Lj in JF IJ — JF so that (24) holds. We first find a basis of C by taking 
repeated Lie brackets of elements in JF. More precisely, set 

:= JF, (39) 

a linearly independent set of elements of 'depth' (no Lie bracket necessary), and let 

Pi:=[Po,^], (40) 

a set of elements of depth 1. which are Lie brackets of elements of depth with elements 
of JF. From the set Vi we extract a possibly smaller set T>i such that Vq U "Z^i is a maximal 
linearly independent set in PqU^'i- Proceeding this way, we now calculate a set of Lie 
brackets of depth 2 

V2:^[Vi,T], (41) 

and extract a subset 'D2 Q 2^2 so that X>o U ^'i U is a maximal linearly independent 
set in T>o[jT>i\JT>2. Proceeding this way, we obtain a set Ufe=o^fe) which spans all of 
C. As a consequence of C being finite dimensional, the procedure will end at some finite 
depth r after which we cannot find any new linearly independent matrix. We write, for 
A; = 0, ...,r, 

'I^k',= {Dlk, D2k, ■ ■ ■ , Dnkk}- (42) 

We can decompose H as 

H = j2Hk, (43) 

k=0 

with Hk a linear combination of elements of depth k, that is, 

nk 

Hk-.^^akjDjk- (44) 
Now, following the proof of the theorem, we can write 

ne^^^ + 0(xi+^). (45) 

k=0 
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Then we can write each of the e as 

^Hkx ^ -Q ^DjkfAx) ^ oix^^^"), (46) 

J=l 

for some Sk > 0. This is straightforward for A; = and it has to be done iteratively for 
Lie brackets of higher depth following the procedure indicated in the proof of theorem. 
Summarizing the method is as follows: 

1. Find a basis for C by repeated Lie brackets of elements of J^. Let r denote the 
maximum depth. 

2. Expand if as a sum of linear combinations of matrices of depth 0, 1, . . ., as in (43), 
(44). 

3. For each of these linear combinations approximate the exponential with a product 
of exponentials involving elements in the basis according to the proof of theorem 
2. In particular the rules to obtain the approximating products are as follows (see 
proof of theorem 2) . 

(a) li A e J^U —T^ then the associated product is T{x) — e^^ (only one factor). 

(b) If T{x) is the product associated with ^4, then T~^{x) is the product associated 
with —A. 

(c) If T{x) is the product associated with A, then T{ax) is the product associated 
with aA for any a > 0. 

(d) If Ta{x) and Tb{x) are the products associated with A and B respectively, 
then Ta{x)Tb{x) is the product associated with A + B. 

(e) If Ta{x) and Tb{x) are the products associated with A and B respectively, 
then TX^{^/x)TB^(^/x)TA{^/x)TB{^/x) is the product associated with 

4. Combine all the products in a unique product approximating e^^, which contains 
only exponentials of elements in JF and — JF. By repeating this product for x — ^ a, 
large number of times n we obtain a matrix arbitrarily close to . 

5. Replace every exponential e"^* with A ^ and t < in the approximating prod- 
uct with an approximating exponential of the form e"^* with i > 0, according to 
proposition 2.1 and remark 2.2. 

3.1 Example 

We illustrate the previous procedure with an example taken from the theory of electrical 
networks. In particular, we consider the LC switching network in [20] (see also [8]) whose 
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dynamical equation is given by 

/O -u 
u 



X 





Vo 










p 
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-p 
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x + 
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V-7 
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xu{t), 



(47) 



where z/, /3, 7 and 5 are positive parameters depending the inductances and capacitances 
of the electrical network. The vector x represents voltages and currents in the network 
and w is a switching control variable which takes values in {0, 1}. To make the discussion 
concrete, we choose the parameters u — 1, P — 3, ^ — 1 and S — 2, so that the set of 
available matrices is 
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The solution of (47) is 

x{t) =X(t)x(0), 
where X — X{t) is the solution of the matrix equation 

X = A{u)X, X(0) = 1, A{1) = Ai, 



A{0) = A2. 



(48) 

(49) 

(50) 



Let us use the notation Ej^ for the skew-symmetric 4x4 matrix which has all the 
entries equal to zero except for the (j7c)-th and (/cj)-th (1 < j < A; < 4) which are equal 
to 1 and —1, respectively. Therefore, we can write 



Ao = -E 
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3E' 



Ai — —E12 + Eu + 2E23 — 3£^34, 
By calculating Lie brackets, at depth 1, we obtain 

A3 := [A2, A,] = -5Ei2 + 7^24, 

at depth 2 

Ai := [A3, Ai] = 17£;i2 + 22£;i4 + 26£;23 + 19^34, and A5 := [A3, A2 



34- 



(51) 



(52) 



22£;i4 + 26^23. 

(53) 



(54) 



At depth 3, we obtain 

Ae := [Ai,Ai] = 145£;i3 - 155£;24. 

As the matrices {Ai}, I = 1, . . . , 6, are linearly independent, they span all of so(4) and 
system (50) varies on the Lie group SO{A), a compact Lie group. 



Let us denote by Tj = Tj{x) the products approximating e j = 1, . . . , 6, and let 



us assume that the control problem is to transfer the state [0, 0, 0, 1]^ to [1, 0, 0, 0]"^. 
choose to drive the transition matrix X in (50) to the value 

/ 
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sin(^) cos(^^'^ 
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(55) 
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We proceed using the composition rules illustrated in (a)-(e) above. Since — [A3, A2], 
we have 

n{x) = Ti\^)T,-'{^)T3{V^)T,{^). (56) 
Moreover, since A3 — [^2,^1] we have 

nix) = T^\^)T{\^)T2{^)T,{^), (57) 

and replacing into (56), we obtain 

n{x) = (58) 

Tfi(^)T2-i(^)Ti(^)T2(^)T2-i(v^)T2-i(^)Tr^(^)T2(^)Ti(^)T2(v/^). 

The product approximating e'^^^i* is T^[^t) which we can express in terms of exponentials 
of Ax and A2 only by replacing Ti and T2 (and and T2^) according to the rules in 
(a) and (b) above. In conclusion, we have 



Q-Ai{^t)^ ^-A2{^t)k ^A^i^t)^ ^A2{^t)i ^-A2{^t)h ^ 



(59) 



We numerically calculated the error 



Err'^in) 



4c -2- 



^ 44 n 



= 8 - 2Tr 



T (-- 
^ U4n 



(60) 



for various values of n and the behavior of the Error as a function of the number of 
iterations n is reported in Table 1. The error goes to zero as predicted by the above 
treatment. In a log-log scale the behavior is essentially linear. 

To conclude the example we have to solve the problem that negative times are not 
allowed and therefore we have to replace terms of the form e'"^^^ and e~^^^, with a; > 
in the expression of T5 with approximations of the form e^'^^ and e"^^^, respectively. In 
the case of A2, since {e'^^*\t e R } is periodic we can always find Xi > such that 
Q-A2X _ gA2a;i every x, and we can simply replace the exponential with x with the 
exponential with Xi in T5, without changing the error. For the exponentials of Ai however 
we need to find an approximation and this is always possible with arbitrary accuracy 
according to proposition 2.1 and remark 2.2. 

To be concrete let us assume that the maximum error we can tolerate is 0.4. From 
Table 1, we choose n = 10^ Fix x ^10"^ We have (cf. (60) and Table 1) 



£;rr( 10^) = e^5 3i _ {T^{x)) 



< 0.235. 



(61) 



Let Ts be the approximation of T5 in (59) where we only use positive values in the 
exponentials, appropriately replacing the exponentials of Ai. In particular, by rewriting 
T^ix) in (59) as 

T,{x) = e-^'^^U,{x)e-^^^^U2{x), (62) 
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Table 1: Results of numerical experiments for the method in section 3. 



number of iterations n 


Error Err 


2 


3.1531 


10 


2.3964 


20 


2.0500 


30 


1.8604 


100 


1.3761 


500 


0.9089 


1000 


0.7599 


5000 


0.5022 


50000 


0.2791 


100000 


0.2341 


500000 


0.1558 


1000000 


0.1301 


5000000 


0.0873 


10000000 


0.0733 


50000000 


0.0490 


100000000 


0.0411 



with Ui{x) = 6-^20:3 gAir.^gA2x3g-A2x2g-A2x3 ^ ^A^x^ ^A^x^ ^A^x^ ^ J^^^g 

n := f,{xi,x) = e^i"ini(a;)e^i"in2(a;). (63) 
Therefore the actual error Err is given by 



Err 



10" 



< 



< 0.235 + 



+ 



[T,{x)r' -[T,{x,x,] 

10^ 



105 



where we used (61). Using the formula for A and B unitary matrices^ 

<n\\A-B\\, 

we write 

Err^< 0.235 + 10^ ||T5(a;) -f5{x,xi)\ . 
^This formula is proved by writing A" — = J2k=i ^"^'^(^ — B)B'^~^, which gives 

n 

m"__B"|| < ^ =n\\A-B\\, 
fe=i 



(64) 



(65) 
(66) 



since multiphcation (right or left) by a unitary matrix does not modify the Probenius norm. 
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In view of our bound on the error of 0.4, we need to find xi > so that T5 (x) — T5 {x, xi ) 
0.165 X 10~^. Now, we have 



< 



Therefore we have 



< 



(67) 
(68) 



+ 



_ ^Ai{x^+xi) 



Therefore, we need to find xi > so that 



|1 _ g^l{a:^+a:i)|| < 



0.165 X 10-^ 



(69) 



We calculate explicitly the eigenvalues of Ai which are given by ±ir and with r and 
I given by 



'l5 + yi25 



'15 - 



(70) 



We have 



||1 - e^i(^^+^i)|| =2^1- cos(r(x3 + xi)) + 1 - cos{l{x^ + xi)). (71) 
Therefore, setting t :— x^ + Xi, formula (69) is certainly satisfied if 

1 - cos(rt) < 8 X 10-^^, (72) 

and 

1 - cos(/t) < 8 X 10"^^. (73) 

Setting e := arccos(l — 8 x 10~^^), we need to find t > x^, positive integers p and q such 
that 

\rt - 2np\ < e, \lt - 27rg| < e. (74) 
One way to do this is as follows. Fix an integer A; > large enough so that 



1 e 



(75) 



According to Dirichlet's approximation theorem (see, e.g., [19] Theorem 1.3) we can find 
p and q, with 1 < p < k positive integers so that 



/ 

-p-q 
r 



1 

^k- 



(76) 
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Choose p and q this way and 

t = ^. (77) 
r 

Using this value of t, the first one of (74) is verified because the left hand side becomes 
zero. Replacing this value of t in the second one of (74) and using (75) and (76) we obtain 
that the second inequality is satisfied as well. Moreover, since ^ > 1, we have that 

t > ^ « 1.7366 >x^ ^ (ii^*^'^) ' ~ 0.0291. (78) 
This concludes the example. 



4 Combination of the two methods 



The main ideas in the two methods of control described in the previous sections can 
be combined in a third method. The main idea of the method in Section 2 was to use 
similarity transformation to generate a basis of the dynamical Lie algebra JC starting from 
the given matrices in in (4) (cf. (9)). The main idea of the method in section 3 is the 
use of the limit in Lemma 3.2, once (24) holds. This allows us to control to the target, 
by repeating a given sequence of available exponentials, with arbitrary accuracy. We can 
combine the two ideas. We first use similarity transformations to obtain a basis of £, 
Ai, . . . , Adim£- Then, if is the target and H — J^jTi^ '^j-^j^ we use the fact that 



Edim C . 
i=i "jA 



dimi2 



n 



(79) 



along with Lemma 3.2 to approximate with arbitrary accuracy the target state, i.e., 

dim C 



lim 

n— >oo 



3=i 



(80) 



At the end of the process, we replace all the exponentials of the form e^^* with t < with 
approximating exponentials of the form e^^* with i > 0. 

We test this method on the example in subsection 3.1. Given Ai and A2 in (48) we 
calculate 



(81) 



1 ° 


-1 





2 ^ 


1 





1 








-1 





-3 


V-2 
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/ 



Our target is e^^3i in (55). We have the decomposition 

^5 = 10^1 + 6F - 16^, 



(82) 
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so that 

e^^M^^R^{x)+0{x^), (83) 

with 

We have, according to Lemma 3.2, 

(85) 

Table 2 shows the results of numerical experiments with this scheme displaying the error 
Err as a function of the number of iterations. Compared with Table 1, it is clear that 
this method converges much faster. Another advantage is that the all the exponentials 
e"^* with negative t are for A = A2 (cf. (84) and (81)) and the one dimensional subgroup 
associated with A2 is closed. Therefore no further approximation is needed. 



lim 

n^oo 



n 



Table 2: Results of numerical experiments for the method in section 4. 



number of iterations n 


Error Err 


2 


2.2819 


10 


0.4544 


20 


0.2267 


50 


0.0906 


100 


0.0453 


1000 


0.0045 



10000 0.0005 



5 Conclusions 

The methods described in this paper can be seen as a constructive proof of the Lie algebra 
rank condition of Theorem 1. It is expected that the ideas described above can be extended 
and improved by using more sophisticated exponential formulas [15], in many ways. It is 
also expected that it will be possible to obtain estimates of the convergence rate in various 
cases. Our goal here was to propose ideas that, although at an early stage, are very general 
and, in principle, allow us to control every system on a compact Lie group. These systems 
in particular include the important class of closed, finite dimensional, quantum systems 
which are coherently controlled, namely controlled through a change in the Hamiltonian. 

In the future, it will be important to improve the algorithms by minimizing the number 
of switches in the control laws that mainly depends on the number of iterations, in the 
last two sections. In this respect, the algorithm of section 4 is expected to be faster than 
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the algorithm in section 3, as a consequence of the exponent 2, in the 0{x'^) in (79) as 
opposed to the exponent 1 + 5 (with S typically < 1) in (24). If our main concern is 
however the time of implementation, the effect of an increasing number of iterations n in 
(26) is balanced by the - exponents inside the limit. The main problem, in terms of time, 
is the approximation of matrices of the form e^^ with t < with matrices of the form 
e"^*, with i > 0, in the case of non-closed subgroups. In fact, we might have to 'travel' 
for a long time inside the Lie group before we get close enough to the original e"^*. 
In special situations, however, it might be possible to transform A into —A via available 
similarity transformations, or reduce ourselves to a smaller dimensional Lie subgroup 
where the problem is more easily tractable. Nevertheless, it is always possible to find 
such an approximation and therefore the control. Remark 2.2 shows how this problem 
can be reduced to a standard problem of Diophantine approximation in number theory 
for which there exist a vast literature and that can be always solved in principle. 

In conclusion, would like to comment on the assumption of compactness which is 
used in the paper only in two instances. In particular, compactness is used to have a 
surjective exponential map and to be able to approximate an exponential of the form e"^* 
with t negative with an exponential of the same type with t positive. Whenever these 
two properties hold, the methods of this paper can still be applied to more general Lie 
groups. In particular this is the case for finite dimensional closed quantum mechanical 
systems whose dynamical Lie algebra £ is a subalgebra of u{n). 
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